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Microwave-induced dynamic switching of magnetic skyrmion cores in nanodots 
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The nonlinear dynamic behavior of a magnetic skyrmion in circular nanodots was studied nu¬ 
merically by solving the Landau-Lifshitz-Gilbert equation with a classical spin model. We show 
that a skyrmion core reversal can be achieved within nanoseconds using a perpendicular oscillating 
magnetic field. Two symmetric switching processes that correspond to excitations of the breathing 
mode and the mixed mode (combination of the breathing mode and a radial spin-wave mode) are 
identified. For excitation of the breathing mode, the skyrmion core switches through nucleation of 
a new core from a transient uniform state. In the mixed mode, the skyrmion core reverses with the 
help of spins excited both at the edge and core regions. Unlike the magnetic vortex core reversal, 
the excitation of radial spin waves does not dominate the skyrmion core reversal process. 


A rich variety of magnetic configurations can emerge 
in nanoscale magnetic materials. For instance, in a cir¬ 
cular nanodisk, both vortex and skyrmion states may 
occur. 1 4 As a topological excitation, a skyrmion has a 
topological charge Q = =bl, whereas the vortex topo¬ 
logical chargd^ is ±1/2. The reversal of a vort ex co re 
has been intensively studied in the past few years 2 4 6 15 
because of its potential use in data storage devices and 
the importance in fundamental physics. Recently, the 
creation and annihilation of individual skyrmions was 
demonstrated,^^ 0 which illustrates the potential appli¬ 
cation of topological charge in future information-storage 
devices. Moreover, the spin-wave modes o f a sk yrmion 
have been inve stiga ted both experimental^ 21 22 and us¬ 
ing simulationsand three k = 0 skyrmion “optical” 
modes have been identified. It was found that anticlock¬ 
wise and clockwise rotation modes are excited when an 
alternating current (AC) magnetic field is applied in the 
skyrmion plane, while the breathing mode is found if 
the AC field is applied perpendicular to the skyrmion 
planeP^ Very recently, a hysteretic behavior and reversal 
of an isolated skyrmion formed in a nanodisk by a static 
magnetic field have been demonstrated.^ However, the 
influence of a time-dependent magnetic field, such as a 
microwave field, on the skyrmion core reversal is still un¬ 
explored. This is the focus of this letter, in which we 
study the reversal of an individual skyrmion in a circular 
magnetic dot by applying an AC magnetic field perpen¬ 
dicular to the skyrmion plane. 

The Dzyaloshinskii-Moriya interaction (DMIp^^ is 
an antisymmetric interaction which ari ses if inversion 
symmetry in magnetic systems is absent plESl either be¬ 
cause of a non-centrosymmetric crystal structure or due 
to the presence of interfaces.^ Consequently, these chiral 
interactions can be classified as bulk or interfacial, de¬ 
pending on the type of inversion symmetry breaking.^ 
In magnetic thin films, a Bloch-type chiral skyrmion 
can be formed by bulk DMI, while a Neel-type radial 
(“hedgehog”) skyrmion requires the presence of inter¬ 
facial DMlP^ In this study we consider the bulk DMI, 
whic h exi sts in non-centrosymmetric magnets such as 

MnSPIHU anc [ p e Q e j33] 


We consider a classical Heisenberg model on a two- 
dimensional regular square lattice with ferromagnetic 
exchange interaction (represented by J) and uniaxial 
anisotropy (K u ) along the z axisP 3 34 l 35 l Apart from that, 
a time-dependent magnetic field h (t) is applied to the 
system in the positive z direction. Therefore, the total 
Hamiltonian of the system is given by 

H = - J ^ m* • nij + ^ T>ij • [m* x rrp] 

<i,j> <i,3> 

i i 

where m* is the unit vector of a magnetic moment 
pi i = — fryS with S being the atomic spin and y(> 0) the 
gyromagnetic ratio. The DMI vector D^ can be writ¬ 
ten as = Dvij in the case of bulk DMI, where D 
is the DMI constant and is the unit vector between 
pi i and pLj. In this study, dipolar interactions are not 
included, because in systems with small sizes of 10 -100 
nm, dipolar interactions are comparably weak. |23 J^ The 
spin dynamics at lattice site i is governed by the Landau- 
Lifshitz-Gilbert (LLG) equation, 

dmi dmi 

—Q t ~ = - 7 m i x H eff ± ami x - 7 ^-, ( 2 ) 

where a denotes the Gilbert damping and the effective 
field H e ff is computed as H e ff = — The Hamil- 

tonian ([l]) associated with the LLG equation Q can be 
understood as a finite-difference micromagnetic model. 
In this study, we have chosen J = 1 as the energy 
unit? 23 1 All simulations are performed in a circular dot 
with a diameter of 121 lattice sites. We have chosen 
J = 7 = 7 = A = las simulation parameters, and 
therefore, the coefficients for conversion of the external 
field h , time £, and frequency uj to SI units are J/(ftr{S ), 
hS/J and J/(hS ), respectively, as shown in Table [l| The 
Gilbert damping a = 0.02 is selected for most of the sim¬ 
ulations. A DMI constant with a value D = 0.08 is used, 
except for Fig. [lj where D is varied. This value of D 
results in the spiral period A « 2tt Ja/D « 39 nm and 
the skyrmion diameter 23 of approximately 55 nm for the 
lattice parameter a = 0.5 nm. 
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TABLE I. Unit 

conversion table for J — 

1 meV and S = 1 . 

Magnetic field h 

J/(trrS) 

« 8.63 T 

Time t 

hS/J 

« 0.66 ps 

Frequency uj 

j/(hS) 

« 1.52 x 10 3 GHz 



FIG. 1. (a) The phase diagram with topological charge of 
the ground state as a function of D and K u . Possible ground 
states are (b) skyrmion (SkX), (c) uniform, and (d)-(e) dis¬ 
torted states. The diameter of the dot is 121 sites. 
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FIG. 2. Imaginary part of the susceptibility spectrum (xzz) 
obtained after applying a field pulse h z = hosmc(uJot) along 
the z axis to the system with ho = 8 x 10 -6 and ujo — 0.1. 
The insets show the spatial distribution of the FFT power at 
eigenfrequencies u — 3.6 x 10 -3 and 1.4 x 10 -2 , respectively. 
The lines are profiles of the m z component across the core 
center. 


There are various stable states for the system de¬ 
scribed by Hamiltonian 0- The phase diagram pre¬ 
sented in Fig. lja) shows the ground state as a func¬ 
tion of D and K u for the magnetic nanodot. The 
phase diagram is obtained by comparing several possible 
states such as skyrmion, uniform, and distorted states, 
as shown in Figs. 0b) -(e). Figure 0b) shows an iso¬ 
lated skyrmion with topological charge (skyrmion num¬ 
ber) Q = f qdxdy « —0.86 where q = m • ( d x m x d y m ) 
is the topological charge density. The isolated skyrmion 
state emerges in the brown region with skyrmion number 
< —0.6, while the distorted states are shown in the purple 
area, where the skyrmion number is chosen to be posi¬ 
tive to distinguish from the skyrmion state. The uniform 
state with skyrmion number « 0 is found in the remain¬ 
ing area. We fix parameters D = 0.08 and K u = 0.004 
in the remaining of this work, as indicated in Fig. []Ja). 

To extract the dynamic response of the skyrmion to 
the microwaves, we calcu late the magnetic absorption 
spectrum of the skyrmion 3 37 by applying a magnetic 
field pulse in the vertical (z) direction with a temporal 
shape given by a sine-function field h z = ho sinc(u;o^) = 
ho sin(cL7 0 t)/(a;ot), where ho = 8 x 10 -6 and ujo = 0.1. 
The magnetic spectrum of a single skyrmion is shown in 
Fig. [2] with a = 0.04. There are two main resonance lines 
at u) = 3.6 x 10 -3 and u) = 1.4 x 10 —2 , corresponding to 
/ = 0.87 GHz and / = 3.39 GHz, respectively, if J = 1 
meV and S = 1 are chosen. The insets in Fig. [2] show 
the corresponding eigenmode amplitude plots, which are 
obtained by computing the fast Fourier transform (FFT) 
to the spatial m z oscillations of the sample. These pat¬ 
terns indicate the spatially resolved fluctuation ampli¬ 
tude. The lower-frequency resonance peak corresponds 
to the breathing mode of the skyrmion, 22 2 ^ where the 


skyrmion radially expands and shrinks as a function of 
time with significant fluctuation around the skyrmion 
core region. The higher-frequency resonance peak with 
weaker intensity can be viewed as a mixed mode that 
combines the breathing mode and the radial spin-wave 
modepl since the spin excitation arises both around the 
core area as well as at the dot edges. 

After determining the eigenfrequencies of the skyrmion 
formed in a nanodot, we study the switching driven by 
a sinusoidal magnetic field along the z direction, i.e., 
h% = h u sin(c ut), where h u and uj are the field amplitude 
and frequency, respectively. Two cases of microwave- 
induced dynamics of a single skyrmion at h u = 2.4 x 10 -3 
(« 20 mT for J = 1 meV and S = 1) with AC fields 
of uo = 3.6 x 10 -3 and 1.4 x 10 -2 corresponding to the 
breathing and mixed resonant modes, respectively, are 
shown in Fig. [3] Figures [3ja) and (b) present serial snap¬ 
shots, while the corresponding m z (t) and h u (t) are shown 
in Fig. |3jc) by black and red curves for the two frequen¬ 
cies. The initial skyrmion state (t = 0) is obtained by 
relaxing the system with maximum \dm/dt\ < 10 -6 . In 
the initial state, the skyrmion number —0.86 and the z 
component of the magnetization m z of the skyrmion core 
are negative, while m z at the edge of the nanodot is pos¬ 
itive. In the first scenario, shown in Fig. [3ja) [See also 
multimedia view of Fig. 3(a)], the breathing mode is 
excited by applying an AC field with uj a = 3.6 x 10 -3 , 
which is switched off at = 1750, as indicated by the 
black dotted line in Fig. |3jc). As a response to the AC 
magnetic field, the skyrmion starts to breathe, i.e., the 
skyrmion core expands when the field points to the neg¬ 
ative ^-direction and shrinks if the direction of applied 
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FIG. 3. Spin dynamics of the skyrmion excited by an AC magnetic held /A = h u ain(ujt) with h u = 2.4 x 10 -3 up to time 
25000 16.5 ns for J — 1 meV and 5 = 1). Selected time evolution of the magnetization distribution with two examples: 

(a) u a = 3.6 X 10 -3 , = 1750 (multimedia view), (b) Ub — 1.4 x 10 -2 , t h 0 ff = 5600 (multimedia view). The color-coded 

m z (out-of-plane) and m x (in-plane) components of the magnetization are displayed in the first and third line, respectively, 
while the second line displays a perspective view (topography) with the color code of m z . (c) Time trace of the magnetization 
component m z (top graph), the AC magnetic field (dashed lines), and the skyrmion number (both in bottom graph). Black and 
red lines refer to the data of the two different examples, while the black circles and red prisms indicate the times corresponding 
to the snapshots shown in (a) and (b), respectively. 


magnetic field is positive. The core thus becomes smaller 
than in the initial state (t = 0) after the first positive half 
period (t = 780). After one full period (t = 1760), the 
skyrmion core expands over the entire dot, the edges be¬ 
come negative, and a transient uniform state is reached. 
Afterwards, the tilted spins at the edge of the disk circu¬ 
late to the inside of the dot and propagate to the center as 
shown in Fig.|3ja) at times 1920-2120. At the same time, 
m z of the core begins to increase and an upward core is 
nucleated at t = 2340 [Fig. |3ja)] . Subsequently, high- 
frequency radial spin waves are excited from the core 
center. The switching of the skyrmion from a downward 
to an upward core orientation is finished. At last the sys¬ 
tem enters into a damped oscillation with the frequency 
of the breathing mode up to t ~ 25000, as shown by the 
black line in Fig. |3jc) . 

Now we turn to the simulation results of the mixed 
mode case, = 1.4 x 10 -2 , with = 5600, shown in 


Fig. ^b) [See also multimedia view of Fig. 3(b)] and by 
the red curves in (c). By applying the AC magnetic field, 
both the spins at the edge of the structure and around the 
skyrmion core are excited in the mixed mode, as shown at 
t = 2720. After t = 3000, the lower-frequency breathing 
mode starts and then dominates the dynamic properties. 
The most expanded core is obtained at t — 3920 with 
a ring of tilted spins around the middle, similar to the 
case of t = 1920 and t = 2120 in Fig. [3](a), however, no 
reversed core is formed for these field parameters. After 
the core has shrunk to a smaller size at t = 4760, it does 
not recover to the previous size because of a compression 
resulting from the excitation of the edge area (t = 5520). 
Instead, the core shrinks even further, and immediately 
thereafter, the down-core vanishes and is replaced by an 
up-core, while spin waves are emitted from the core area 
(see t = 5720). After switching off the AC field at t = 
5600, the reversed skyrmion starts to dissipate the energy 
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via both the breathing and mixed modes. 

Si milar to t he vortex core reversal with out-of-plane 
fieldsJ^ 3 15 ^ 3 ^the reversal processes for the two frequen¬ 
cies are symmetric, and the switching mechanisms are 
similar to the radial-spin-wave-mode-assisted switching 
for a vortex core! 13 At the lower frequency, which corre¬ 
sponds to the breathing mode, the skyrmion first switches 
to a transient uniform state via axial core expansion to 
the edge, then a new core is formed from the center. For 
the mixed mode excitation at the higher frequency, the 
reversed skyrmion results from the shrinking of the core 
with the help of spin excitations from the edge region. 
High-frequency spin waves are emitted during the core 
reversal, however, a spin-wave assisted re-reversal, i.e., a 
switching back to the original state by a spin wave, as it 
has been observed in vortex core reversal)- 13 15 38 has not 
been found here. Our simulation results also show that 
the switching time and switching field for a skyrmion- 
core reversal are similar to those needed for a vortex-core 
reversal. 
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FIG. 4. Phase diagram of the skyrmion switching induced by 
a field h% — h^sin^t) as a function of h V and u. Red dots 
indicate skyrmion reversal, gray dots mean no switching, up- 
pointing and down-pointing triangles skyrmion switching to 
uniform and distorted states, respectively. 


Besides the duration time, the skyrmion reversal pro¬ 
cess also depends on the frequency uj and the AC field 
amplitude h u . Figure [4] shows the switching phase dia¬ 
gram as a function of uj and h u with a maximum temporal 


duration t = 10000 6.6 ns for J = 1 meV and S = 1). 

The diagram can be divided into three different regions. 
The first is the weak-susceptibility region I (Gray sym¬ 
bols), in which a skyrmion core never switches for AC 
field duration up to t = 10000, mainly because the field 
is not strong enough with respect to the dynamic sus¬ 
ceptibility of the skyrmion. In the second region (II) 
a skyrmion converts to either the uniform state (up- 
pointing triangles) or the distorted state (down-pointing 
triangles). The distorted state is formed directly from the 
skyrmion due to the instability of the core under the AC 
field with intermediate intensities, similar to the case of 
magnetic vortices The boundary between regions I and 
II reflects the susceptibility spectrum of the skyrmion, cf. 
Fig-0 in such a way that at around the eigenfrequencies 
of the system switching to the uniform or distorted states 
occurs at lower field amplitudes than at other frequen¬ 
cies. In the third region (III, red), a skyrmion core can 
be reversed by the AC field, as had been shown for two 
instances with = 1750 and t h oE = 5600 in Fig. [3](a) 
and Fig. !b), respectively. Comparing to the results 
for a static magnetic field (uj = 0), where in this ampli¬ 
tude regime the final state is either unchanged or uni¬ 
form depending on whether is below or above about 
3.6 x 10 -3 , it is obvious that an AC magnetic field within 
a certain frequency range substantially helps to reverse 
the skyrmion core. 

In summary, we have studied the nonlinear skyrmion 
dynamics in a circular magnetic nanodot induced by a 
vertical oscillating magnetic field using micromagnetic 
simulations. We found that in a certain frequency win¬ 
dow a fast skyrmion switching can be achieved for field 
amplitudes that do not lead to a reversed skyrmion un¬ 
der static conditions. We presented two examples of 
skyrmion switching in detail, in which the excitation fre¬ 
quencies corresponded to the breathing and mixed modes 
of the system. Under excitation with the breathing-mode 
frequency, the skyrmion reverses via a transient uniform 
state, while exciting it with the mixed-mode frequency 
spin waves from the edge area assist the switching. Our 
results show that a skyrmion core can be reversed within 
nanoseconds by means of microwaves perpendicular to 
the skyrmion plane. 
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